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Abstract 

Single  surface  multipactor  in  high  power  microwave  systems  was  investigated  com¬ 
putationally  and  analytically.  The  research  focused  upon  understanding  the  cause  and 
parametric  dependence  of  the  multipactor  process  leading  to  suggested  methods  of 
mitigation.  System  damage  due  to  reflection  was  also  assessed.  All  simulations  were 
performed  using  the  PIC  code  developed  by  AFRL,  known  as  ICEPIC.  In  recreat¬ 
ing  the  susceptibility  curves  that  define  regions  of  multipactor  growth  and  decay,  a 
discrepancy  was  found  between  previous  published  results  and  those  observed  in  the 
current  simulation.  This  was  attributed  to  previous  simulations  not  accounting  for 
the  magnetic  component  in  the  electromagnetic  radiation  incident  on  the  dielectric 
window.  By  surveying  different  static  magnetic  and  electric  fields  both  parallel  and 
perpendicular  to  the  dielectric,  revised  susceptibility  curves  were  determined.  An 
analytic  method  confirmed  the  origin  of  the  discrepancy.  The  reflection  of  the  waveg- 
uided  radiation  by  the  multipactor  electrons  degrades  the  output  and  may  damage  the 
microwave  source.  A  theory  for  the  reflection  resulting  from  multipactor  was  derived 
to  aide  in  quantifying  reflection.  A  plane  wave  was  used  to  model  incoming  microwave 
radiation  and  the  amplitude  of  the  electric  field  in  reflected  and  transmitted  waves 
was  measured  for  various  frequencies  of  the  incident  plane  wave. 
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MULTIPACTOR  DISCHARGE  IN  HIGH  POWER  MICROWAVE  SYSTEMS: 
ANALYZING  EFFECTS  AND  MITIGATION  THROUGH  SIMULATION  IN 

ICEPIC 


I.  Introduction 


1.1  Purpose 

The  DoD  interest  in  researching  microwave  technology  began  in  the  early  1930s 
and  substantially  increased  during  World  War  II  due  to  the  rise  of  radar  technology 
utilized  by  the  Allied  powers  throughout  the  war  [4],  After  the  war,  many  of  the 
processes  developed  with  higher  power  radar  systems  led  to  many  commercial  uses  of 
microwaves  in  fields  of  communication,  navigation,  spectroscopy  and  radiative  power 
transfer  [18]  [7].  Even  though  the  private  sector  has  made  ample  use  of  microwave 
technology  since  the  end  of  World  War  II,  the  Department  of  Defense  still  employs 
research  in  the  development  of  microwave  systems  in  the  furtherance  of  tactical  and 
strategic  capabilities  [4,  p.  8].  In  1975,  experiments  conducted  by  both  the  Naval 
Research  Laboratory  and  Cornell  led  to  the  development  of  free  electron  MASERS 
capable  of  generating  1  GW  in  the  X-band  of  the  microwave  spectrum 

The  most  recent  DoD  research  regarding  microwave  weapons  follows  along  the 
lines  of  the  Air  Force  Scientific  Advisory  Board’s  New  World  Vistas  statement: 


“Promising  present-day  research  in  high  power  microwave  (HPM)  tech¬ 
nology  allows  us  to  envision  a  whole  new  range  of  compact  weapons  that 
will  be  highly  effective  in  the  sophisticated,  electronic  battlefield  environ¬ 
ment  of  the  future.  There  are  many  advantages  of  these  HPM  weapons: 
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First,  there  are  virtually  no  target  acquisition,  pointing,  and  tracking  re¬ 
quirements  in  any  HPM  weapon  employment  scenario.  Electromagnetic 
(EM)  radiation,  traveling  at  the  speed  of  light,  will  envelop  a  large  volume 
that  can  engage  multiple  targets  at  once.  HPM  weapons  could  be  used  in 
nearly  all  weather,  although  frequencies  above  10  GHz  degrade  somewhat. 
Some  HPM  weapons  designs  could  be  employed  in  a  covert  way  since  the 
beam  is  not  visible  and  the  damage  and/or  upset  could  be  directed  to 
electronic  targets.  Since  in  many  applications  the  only  expendable  is  fuel 
for  electrical  generators,  HPM  weapons  are  expected  to  come  with  a  large 
“magazine.”  It  may  be  possible  to  design  a  system  that  acts  as  both  radar 
and  weapon,  which  Erst  detects  and  tracks  the  target  than  increases  the 
power  and  engages  the  target,  all  at  electronic  speeds,  Finally,  a  distinct 
advantage  of  HPM  lies  in  the  fact  that  it  may  be  considered  a  nonlethal 
weapon  that  would  prevent  the  enemy  from  using  his  electronic  equipment 
successfully  with  no  impact  on  human  life.  [1].” 


This  quote  illustrates  the  importance  of  HPMs  to  modern  military  operations. 
However,  it  does  not  mention  the  technological  hurdles  that  must  first  be  crossed  be¬ 
fore  such  benefits  can  be  realized  by  the  modern  warfighter.  The  Defense  Technology 
Area  Plan  lists  these  required  developments  for  the  success  of  HPM  weapons  [2]: 

1.  Compact,  high-peak-power  or  high-average-power  HPM  sources. 

2.  Compact,  high-gain,  ultra-wideband  antennas 

3.  Compact,  efficient,  high-power, pulsed  power  drivers. 

4.  Compact,  efficient  prime  power  sources. 

5.  Predictive  models  for  HPM  effects  and  lethality 

6.  Low-impact  hardening  of  systems  hostile  and  self-induced  Electro-Magnetic 
Interference 

7.  Affordable  system  integration  meeting  military  platform  requirements. 
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1.2  Problem  Description 


This  research  focused  on  improving  items  1  and  2  of  the  technology  plan  by 
analyzing  the  source  initiated  secondary  electron  emission  from  the  dielectric  window 
that  separates  the  interior  of  the  HPM  waveguide  from  the  atmosphere.  Thermal 
and  mechanical  window  failure  compromises  the  vacuum  integrity  of  the  microwave 
system  and  limits  the  maximum  output  power. 

Current  thought  concerning  HPM  window  damage  indicates  the  cause  stems  from 
thermodynamic  effects  due  to  the  formation  of  multipactor  along  the  cell  window  [13] . 
Multipactor  is  characterized  by  the  creation  of  an  electron  cascade  along  the  inside  of 
the  waveguide  window  in  the  presence  of  high  amplitude  microwave  radiation,  thereby 
creating  plasma  inside  the  waveguide  vacuum.  The  presence  of  plasma  not  only 
damages  the  window  through  thermodynamic  effects,  but  also  leads  to  second  order 
effects  that  end  up  damaging  the  entire  HPM  system.  For  example,  the  reflection  of 
EM  radiation  by  the  multipactor  tends  to  damage  the  HPM  source  and  lowers  the 
power  output  of  the  entire  system. 

1.3  Goals  of  the  Research 

This  research  focused  on  the  characterization  of  multipactor  with  the  intended 
end  result  of  identifying  methods  of  mitigation  of  secondary  electron  growth.  In 
order  to  work  towards  these  end  results,  multipactor  initiation  was  surveyed  to  better 
understand  under  what  circumstances  multipactor  growth  and  decay  occur.  Power 
balance  was  analyzed  and  an  analytic  development  of  reflection  was  developed 
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1.4  Background 


In  1934  Otto  Farnsworth  was  the  first  to  recognize  the  phenomenon  of  multipactor 
discharge.  He  named  the  phenomenon  after  the  effects  he  observed  in  his  invention, 
the  “AC  electron  multipactor”  [19].  Multipacting  was  also  observed  in  many  other 
gaseous  breakdown  experiments  and  was  characterized  by  Abraham,  Sperry  Gyro¬ 
scope  engineers  ,  and  by  Forrer  and  Milazzo  in  1950’s  [23].  Due  to  the  lack  of  in¬ 
formation  regarding  the  dynamics  of  secondary  electron  emission,  Vaughan  extended 
theoretical  knowledge  of  secondary  electron  emission  with  empirical  formula  charac¬ 
terizing  electron  emission  from  a  surface  as  a  function  of  impact  electron  energy  and 
angle  [24], 

With  Vaughan’s  empirical  model,  the  process  behind  multipactor  becomes  clear. 
The  multipactor  process  follows  from  an  initial  seed  electron  impacting  the  dielec¬ 
tric  and  causing  secondary  electron  emission.  These  secondary  electrons,  then  being 
ejected  from  the  surface,  are  accelerated  by  the  parallel  electric  field  associated  with 
microwave  radiation.  The  electrons  are  then  returned  to  the  surface  by  either  a  per¬ 
pendicular  bias  electric  field,  the  positive  charge  on  the  dielectric  window  left  by  the 
secondary  emitted  electrons,  or  the  magnetic  field  component  of  the  microwave  radia¬ 
tion  propagating  through  the  window.  The  acceleration  from  the  microwave  radiation 
results  in  an  increase  in  the  electron’s  kinetic  energy  during  it’s  flight  in  the  waveg¬ 
uide  before  impacting  the  window.  If  the  particle’s  kinetic  energy  is  great  enough, 
the  impacting  electron  will  cause  more  electron’s  to  be  emitted.  For  high  electric 
field  amplitudes  this  process  results  in  an  avalanche  of  electrons  on  the  outside  of  the 
waveguide  dielectric  window.  An  illustration  of  this  process  can  be  seen  in  Figure  1. 

Previous  work  in  this  field  has  been  focused  on  the  analysis  of  window  breakage  in 
HPM  systems,  with  experimental  analysis  performed  at  Texas  Tech  University.  These 
experiments  focused  on  characterizing  dielectric  window  breakdown  in  the  presence 
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Figure  1.  The  multipactor  effect  on  a  waveguide  window.  Electrons  are  emitted  from 
the  vacuum  side  of  the  dielectric  window  and  are  pulled  back  to  the  window  due  to  the 
formation  of  a  positive  charge  (with  a  field  labeled  as  E± )  and  gain  energy  from  the 
electromagnetic  field  (labeled  as  E ||). 

of  high  intensity  microwave  radiation  [4],  The  laboratory  experiments  performed 
by  Texas  Tech  also  observed  the  multipactor  discharge  along  the  dielectric  window 
through  x-ray  emission,  discharge  luminosity,  and  microwave  fields  [16] .  Other  studies 
performed  over  the  past  few  years  have  also  indicated  that  the  magnetic  fields  play  a 
large  effect  in  determining  electron  transit  time  [11]. 

Computational  work  has  also  joined  the  research  effort,  being  cost  effective  and  a 
vehicle  to  understanding  the  fundamental  physics  involved  in  the  multipactor  process. 
Monte  Carlo  techniques  were  originally  used  to  process  the  electron  avalanche  through 
the  simulation  of  an  oscillating  E\\  held  representing  the  electric  held  component  of 
incoming  microwave  radiation,  and  a  static  E±  held  representing  either  the  return 
held  associated  with  space  charge  present  on  the  dielectric  window  or  an  imposed  bias 


5 


[13].  Electron  avalanche  was  accurately  modeled  with  Vaughan  formula  for  secondary 
electron  yield  [8].  Experimental  research  performed  at  Texas  Tech  [16]  showed  that 
space  charge  must  be  included  to  properly  characterize  the  multipactor  evolution  [22] . 
This  has  led  to  the  use  of  Particle-in-Cell  (PIC)  methods  for  modeling  of  dielectric 
window  systems  [19]  [20]. 

1.5  Definition  of  Terms 

Before  delving  into  the  research  accomplished,  it  is  important  to  characterize  two 
of  the  different  situations  observed  throughout  the  simulations.  One  in  which  space 
charge,  the  charge  associated  with  the  electrons  and  the  positive  charge  contribution 
from  the  dielectric,  is  ignored.  This  case  will  be  referred  to  as  the  Non-Space  Charge 
(NSC)  case.  The  Space  Charge  (SC)  case  is  characterized  by  the  allowance  of  particle 
growth  and  decay  through  the  self-consistent  inclusion  of  the  electron’s  effects  on  each 
other  and  the  effects  of  the  positive  charge  on  the  surface  of  the  dielectric. 

Before  proceeding,  basic  definitions  of  relevant  physical  quantities  and  their  cor¬ 
responding  symbols  are  provided. 

•  An  E  symbol  stands  for  the  electric  field  and  B  stands  for  the  magnetic  field. 

•  An  £  indicates  energy  associated  with  the  free  electrons 

•  An  u  indicates  angular  frequency. 

•  The  physical  constants  of  mass  and  charge  are  written  as  m  and  q  respectively. 

•  These  definitions  may  be  further  modified  by  the  use  of  subscripts  indicating 
the  association  of  the  physical  quantity  with  either  a  direction  or  object  within 
the  system.  For  example,  the  symbol  qe  indicates  the  charge  of  an  electron  and 
the  symbol  vz  indicates  the  component  of  thevelocity  in  the  z  direction. 

•  Some  of  the  symbols  used  throughout  the  paper  will  be  different  than  the 
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conventional  ones  used  throughout  much  of  the  literature  due  to  improvements 
in  the  understanding  of  the  multipactor  process. 

Throughout  the  paper  EFiF  held  symbol  will  be  replaced  with  the  symbol 
En  in  representing  the  E  held  parallel  to  the  window.  This  admits  treat¬ 
ments  of  oblique  incidence  where  the  EFiF  held  contributes  components 
parallel  and  perpendicular  to  the  dielectric  window. 

The  Epc  will  be  replaced  with  the  symbol  E±  ,  representing  the  E  held 
perpendicular  to  the  window.  This  change  is  due  to  the  the  electric  held 
associated  with  the  positive  charge  on  the  dielectric  window  dynamically 
evolving  as  a  function  of  particles  in  the  SC  case  and  accommodates  the 
oblique  incidence  condition 

Further  definitions  of  symbols  will  be  discussed  as  they  occur  in  the  paper. 

1.6  PIC  Code 

This  section  outlines  the  use  of  the  PIC  code  used  in  modeling  the  movements 
of  particles  within  the  system  used  throughout  all  of  the  research.  Particlc-in-ccll 
methods  work  by  tracking  individual  particles  in  a  continuous  space,  but  treats  the 
calculations  of  their  held  elements  as  functions  of  both  current  and  particle  densities 
per  cell.  This  method  is  efficient  because  by  representing  a  large  number  of  particles 
as  a  single  large  pseudo  particle,  the  physics  of  a  problem  may  be  better  understood 
while  increasing  the  efficiency  of  a  simulation.  The  typical  PIC  algorithm  operates  in 
a  four  step  cycle,  in  which  the  particle,  held  interactions,  and  energy  exchanges  are 
treated  self-consistently.  An  illustration  of  the  cycle  can  be  found  in  Figure  2.  For 
this  research,  the  Improved  Concurrent  Electromagnetic  Particle-in-Cell  (ICEPIC) 
was  chosen  in  modeling  the  system.  ICEPIC  is  a  self  consistent  and  fully  relativistic 
modeling  software  developed  by  Air  Force  Research  Laboratory  (AFRL)  [5]  to  aid  in 
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Figure  2.  The  four  step  PIC  algorithm  illustrating  the  process  of  taken  by  the  PIC 
code  for  every  iteration. 


the  development  of  HPM  devices. 


1.7  Consistency  of  PIC  Methods 

Even  though  PIC  methods  grant  insight  of  the  physics  behind  a  system,  PIC 
operations  require  that  a  few  conditions  must  be  met  in  order  for  accurate  results 
of  systems  to  be  observed.  One  such  constraint  involves  the  time  step  (dt)  of  the 
simulation.  Normally,  the  thermal  velocity  of  the  particles  necessitates  the  size  of 
the  time  step  must  be  less  than  where  dz  is  the  spatial  size  of  the  grid  and  v 
is  the  velocity  of  thermal  velocity  of  an  electron.  When  considering  fully  coupled 
electrodynamic  simulations,  the  requisite  velocity  is  the  speed  of  light.  Therefore,  the 


simulation  timestep  is 


(1) 


dt  = 


0.99  dz 


where  c  is  the  speed  of  light. 

Another  constraint  requires  that  the  grid  scale  resolve  structures  smaller  than  the 
Debye  length  associated  with  the  multipactor  plasma  and  is  complemented  by  the 
additional  statistical  requirement  that  each  cell  in  the  simulation  containing  at  least 
5-10  particles  [6].  Previous  results  in  meeting  these  requirements  have  yielded  data 
that  agree  with  analytic  theory  [13]  [19]. 


Simulation  Setup. 

The  default  material  of  the  ICEPIC  environment  is  a  perfect  conductor.  Vacuum 
within  the  system  must  be  carved  out  of  this  conductor  allowing  for  the  simulation  of 
an  environment  comparable  to  the  inside  of  a  microwave  waveguide.  The  simulated 
geometry  model  used  follows  the  model  described  both  by  Fichtl  [8]  and  Rogers  [19] 
and  is  displayed  in  Figure  3  . 

The  Perfectly  Matched  Layer  (PML)  of  the  simulation  acts  as  perfect  absorbers 
of  the  plane  wave  moving  through  the  waveguide.  Another  important  feature  in  the 
simulation  is  the  inclusion  of  a  secondary  dielectric  that  prevents  emitted  high  energy 
particles  from  interacting  with  the  plane  wave  emitter  on  the  —z  side  of  the  simulation. 
Finally,  a  periodic  boundary  condition  is  also  included  along  the  waveguide,  meaning 
that  any  particles  traversing  the  boundary  shall  return  on  the  other  side  of  the  periodic 
boundary  of  the  system.  This  allows  for  a  smaller  simulation  grid-space  to  be  used, 
while  still  maintaining  the  particle  rich  environment  of  the  simulation. 

In  order  for  the  system  to  maintain  a  consistent  resolution  while  varying  the  fre¬ 
quencies  used  in  the  simulation,  dz  is  made  to  be  dependent  on  wavelength  divided 
by  a  resolution  constant  (factors  of  3000  for  when  space  charge  is  not  included  and 
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Figure  3.  Simulation  geometry  following  Fichtl  [8]  and  Rogers  [19]. 


12000  for  when  space  charge  is  included  in  the  simulation).  This  difference  in  resolu¬ 
tion  was  based  upon  previous  work  done  by  Fichtl  [8],  and  comes  from  the  previous 
requirement  to  maintain  the  particle  dynamics  by  constraining  the  grid  size  to  be 
smaller  than  a  Debye  length  [8]. 

For  most  ICEPIC  simulations,  the  frequency  of  the  sinusoidal  oscillations  of  the 
electric  and  magnetic  fields  was  2.45  GHZ .  For  the  case  where  space  charge  is  treated, 
nominal  specifications  of  the  system  were 
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dz  = 

A  ft*  1.02  x  10“5m 

(2) 

12000 

A.X  Full  ~ 

Nx  x  dz  ft*  7.10  x  10_5m 

(3) 

Az Particle  Space 

156  x  dz  ft*  1.59  x  10_3m 

(4) 

AzpML  Depth  = 

150  x  dz  ft*  1.53  x  10_3m 

(5) 

Az Simulation  Space 

Nz  x  dz  ft*  5.10  x  10_3m 

(6) 

as  indicated  in  Figure  3. 

The  temporal  resolution  of  the  system  depends  on  the  value  dz,  so  that  the  tem¬ 
poral  resolution  may  be  treated  as  a  function  of  frequency  and  resolution, 


dt 

dt 

dt 


0.99  dz 

V2c 

0.99A 

12000^  c 

0.99  14 

- «  2.38  x  10~14s 

12000  y/2f 


(7) 

(8) 
(9) 


This  time  step  value  is  roughly  a  factor  of  17000  smaller  than  the  period  of  the 
incoming  wave,  and  on  the  order  of  one  hundredth  the  transit  time  of  the  electron. 
Therefore,  the  simulation  accurately  model  the  physics  associated  with  a  real  mi¬ 
crowave  system. 
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II.  Understanding  Multipactor:  Susceptibility  Mapping 


A  susceptibility  diagram  is  a  parametric  map  of  regions  of  growth  and  decay 
in  the  space  defined  by  E\\  and  Ej_ .  This  is  a  logical  choice  because  the  injected 
electron  gains  energy  primarily  from  the  parallel  component  of  the  microwave  field 
while  the  perpendicular  component  determines  the  time  of  flight.  In  mapping  the 
susceptibility  diagram  for  multipactor  initiation,  the  simulations  focused  on  the  NSC 
case.  Considering  the  differences  in  implementation  of  the  NSC  and  SC  in  both  theory 
and  in  the  ICEPIC  simulation  it  is  convenient  to  treat  the  mapping  of  the  initiation  of 
multipactor  growth  and  decay  in  this  chapter  and  reserve  the  analysis  of  multipactor 
dynamics  and  power  balance  to  Chapter  III. 


2.1  Theory 


Secondary  Electron  Emission. 

Vaughan’s  emission  curves  allow  for  empirical  modeling  of  the  secondary  electron 
yield  from  data  produced  by  others,  such  as  Gibbons  [24],  The  geometry  of  the 
dielectric  surface  on  which  this  modeling  is  based  is  shown  in  Figure  4. 

For  this  configuration  the  ratio  of  secondary  emitted  electrons  to  incident  electrons 
(  5  )  in  terms  of  the  angle  of  impact  (  <^  )  is 


Smaxk)  ■  {w  exp((l  -  w)))k{-w)  for  w  <  3.6 

<*(ei>s)  =  <  (10) 

•  (1.125-u;-0-35)  for  w  >  3.6. 


Where  w  =  — and  the  corrected  maximum  electron  emission  energy  is  given  by 

£max  ) 


ksS2 

£ maxis )  ^max  0  (1  T  )• 

7 T 


(11) 
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Figure  4.  Trajectory  of  the  secondary  electron,  vo  is  the  electrons  initial  energy,  <j>  is 
the  electron’s  initial  emission  angle,  and  <;  is  the  angle  at  which  the  particle  impacts 
the  dielectric  surface. 

The  corrected  maximum  emission  term  in  respect  to  the  impact  angle  to  the  surface 
may  be  written  as 

k  c2 

dmax(s)  =  Smax o  •  (1  +  — — )•  (12) 

Z7T 

In  both  of  these  equations  ks  is  the  smoothness  parameter  of  the  system  (usually 
set  equal  to  1  [19]).  emaxo  and  5max0  are  constants  dependent  on  properties  of  the 
material  and  may  be  altered  to  effect  secondary  electron  emission.  In  Equation  10, 
k{w)  is  also  normally  defined  as  either  as  a  piecewise  function  of  w: 

f0.56  for  w  <  1 

(13) 

0.25  for  1  <  w  <  3.6 
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Figure  5.  Emission  yield  vs  Impact  energy.  The  various  curves  represent  different 
values  for  5maxo,  an  input  based  on  the  dielectric  window  material.  e1  and  £2  are  also 
known  as  the  low  and  high  crossover  energies.  This  figure  is  based  on  the  figures 
generated  from  Vaughan’s  equations[24]. 


or  in  terms  of  Valfells  continuous  representation  of  k  in  terms  of  w  [22]: 


k  =  0.435 


0.27 
7 r 


arctan(-7T  ln(iu)). 


(14) 


Figure  5  shows  an  example  output  of  the  function  for  different  values  of 

Smax 0,  indicating  the  degree  to  which  the  secondary  electron  yield  depends  upon  the 
material  properties  of  the  dielectric.  The  figure  illustrates  that  the  energy  associated 
with  electron  impact  dictates  the  secondary  emission  yield  associated  with  particle 
impact.  When  5  >  1  an  electron  impacting  on  the  surface  indicates  that  more  than 
one  electron  is  emitted  for  the  impact.  Likewise,  when  5  <  1  the  impacting  electron  is 
absorbed.  Within  the  theoretical  analysis,  it  was  assumed  that  the  impact  angle  is  | 
due  to  the  large  velocity  component  expected  for  motion  in  the  x  direction  due  to  the 
accelerating  E»  held.  This  means  that  for  the  research,  only  two  impact  energies  were 
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Figure  6.  A  colored  secondary  emission  curve,  where  the  color  green  indicates  im¬ 
pact  energies  that  lead  to  secondary  particle  growth  and  the  color  red  corresponds  to 
energies  at  which  secondary  particle  growth  will  not  occur. 

considered  for  analytically  determining  whether  the  the  impacting  electron  generates 
more  electrons  or  gets  absorbed.  Therefore,  the  lower  cutoff  energy  (labeled  ey)  and 
upper  cutoff  energy  (labeled  e 2 )  determine  a  region  of  multipactor  growth  and  decay. 

Initiation  of  Multipactor. 

The  secondary  emission  curve  indicates  boundaries  governed  by  the  average  energy 
of  particle  impact.  If  the  average  energy  of  the  particles  impacting  the  surface  of 
the  dielectric  is  greater  than  the  energy  at  ey  and  below  the  higher  cutoff  energy 
of  £2  multipactor  will  be  initiated  within  the  system.  Otherwise  electrons  will  be 
absorbed  and  multipactor  shall  cease.  These  regions  of  the  secondary  emission  curve 
are  color  coded  in  Figure  6.  With  the  the  microwave  field  (E»  )  acting  as  the  sole 
energy  provider  to  the  electron,  it  becomes  evident  that  multipactor  initiation  may 
be  written  in  terms  of  the  En  field.  However,  the  E±  field  must  also  be  considered 
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for  it’s  role  in  returning  the  particle  to  the  dielectric  window.  By  writing  the  initial 


energy  of  the  electron  as  a  function  of  velocity  (emo  =  ^mevf)  and  considering  only 
the  E±  as  the  return  force(i^  =  qeE±),  the  time  of  flight  [22]  of  the  particle  may  be 
written  as 

V2em0  me 

T  =  - — - 

qeE± 

Assuming  the  energy  gained  is  proportional  to  the  time  of  flight,  the  secondary  elec¬ 
tron  yield  may  be  written  as  a  function  of  both  E\\  and  E±  .  This  redefinition  of 
secondary  electron  yield  as  being  parametrically  defined  by  the  electric  fields  parallel 
and  perpendicular  to  the  dielectric  surface  yields  a  susceptibility  diagram.  An  exam¬ 
ple  susceptibility  diagram  is  shown  in  Figure  7.  Note  that  within  the  diagram  the 
region  of  multipactor  growth  (shown  in  the  figure  as  a  green  region)  corresponds  to 
energies  between  crossover  energies  E\  and  e 2 .  The  red  region  in  the  susceptibility 
diagram  correspond  to  regions  of  particle  decay.  The  lines  on  the  example  figure 
correspond  to  boundary  curves,  these  are  regions  where  the  secondary  electron  yield 
is  equal  to  one.  5  =  1  defines  the  boundary  curves  at  which  multipactor  growth  or 
decay  occur.  These  boundary  lines  are  important  in  determining  the  degree  to  which 
a  bias  held  needs  to  be  utilized  to  prevent  multipactor  initiation. 

Analytic  Solution  for  Determining  Boundary  Curves  with  an  Oscillating 

Electric  Field. 

Previously,  boundary  curves  have  been  created  using  the  analytic  expressions  orig¬ 
inally  formulated  by  Valfell’s  in  his  initial  treatment  of  characterizing  multipactor 
growth  [22] .  This  may  be  accomplished  by  deriving  the  electron  energy  gained  during 
its  flight  in  the  E\\  held.  By  treating  the  E\\  held  as  a  function  E\\  *  sin  (cut  +  6)  where 
co  is  the  angular  frequency  of  the  EM  wave,  6  is  the  emission  phase,  and  E n  is  the 
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Decay 


Ei  (  MVj  ml 

Figure  7.  Susceptibility  diagram  for  the  NSC  case.  The  regions  are  colored  based  on 
whether  the  region  will  initiate  multipactor.  Red  indicating  that  multipactor  shall  not 
be  initiated  and  green  indicating  that  multipactor  shall  be  initiated.  These  regions 
correspond  to  average  impact  energies  graphed  in  Figure  6. 


amplitude  of  the  electric  field,  the  equations  of  motion  for  the  electron  become 


y 

z 


eE\\ 

me 


sin  [ut  +  6] 


me 


(16) 

(U) 


The  energy  gained  by  an  electron  during  a  complete  time  of  flight  is 


Aey  =  -my2(r) 


Oy 


(18) 


where  Aey  is  the  change  in  kinetic  energy  of  the  electron,  Vqv  is  the  electrons  initial 
velocity  perpendicular  to  the  surface  of  the  dielectric,  and  r  is  the  time  of  flight.  By 
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combining  this  with  the  solution  to  equation  16,  Equation  18  can  be  rewritten  as 
i  E  e2E2 

Aev  =  —iA — II  [cosfcar  +  6)  —  cos(d)]  H - L  [cosftar  +  6)  —  cos(d)]2  .  (19) 

uj  2  meuj2 

When  Equation  19  is  averaged  based  on  the  emission  phase,  the  average  energy  be¬ 
comes 

^  I1  -  cos<WT>i  <20> 

where  the  electron  time  of  flight  (r)  was  previously  indicated  in  Equation  15.  As¬ 
suming  mono-energetic  emission,  where  the  electron  emission  energy  equated  to  the 
average  of  the  emission  energy  distribution  [8] ,  the  electron  transit  time  then  becomes 


Ay/nieSom 

T  =  — 1 - - . 

qeE± 

Inserting  this  value  into  Equation  20  yields  the  full  energy  gain 


(A^> 


2  rneu2 


-1 

4  y/me£0m 

1  —  cos 

QeE± 

UJ  j 

(21) 


(22) 


and  by  adding  this  equation  to  the  original  injection  energy  (e0m)  leads  to  an  expres¬ 
sion  for  the  final  impact  energy  on  the  dielectric 


£  Impact,  ^Om  "T 


„2  7712 

qeE\\ 

2meu>2 


1  —  cos 


4  y/me£om 


qeEE , 


(23) 


This  may  be  rewritten  as 


N 


2  'TTle£  Impact  2£q  mX^e 


1  —  COS 


4uiy/mee  o  m 
qeE± 


qeE\ 
L 0 


(24) 


18 


By  setting  £/mparf  to  crossover  energies,  the  boundary  curves  associated  with  5=1 
may  be  determined.  Equation  24  may  be  simplified  by  assuming  yjme£ omu;  <C  qeEoc 
leading  to 

Edc  =  Erf-  (25) 

This  yields  the  linear  relationship  previously  observed  in  boundary  curves  for  the 
NSC  cases  reported  by  Ficthl  and  Rogers  [8]  [19]. 

Analytic  Solution  for  Determining  Boundary  Curves  with  Static  Mag¬ 
netic  and  Electric  Fields. 

In  order  to  corroborate  findings  found  in  remapping  the  susceptibility  diagrams, 
it  became  important  to  also  observe  the  effect  of  a  magnetic  field.  Understanding 
that  the  electron  transit  time  is  much  shorter  than  the  period  of  the  electromagnetic 
radiation,  one  may  treat  all  of  the  fields  associated  with  microwave  energy  as  being 
static  where  the  force  equation  is 

F  =  qeE  +  qe(v  x  B).  (26) 

For  this  analytic  development,  the  electric  and  magnetic  fields  associated  with  the 
electromagnetic  radiation  will  be  written  as 


B\\ 

—  brf 

(27) 

E\\ 

=  Erf 

(28) 

and  the  bias  field  shall  be  written  as 

E± 

=  —Eric- 

(29) 
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The  equations  of  motion  associated  with  the  trajectory  of  the  electron  are 


d2y 
dt 2 
d2z 
dt2 
dy_ 

dt 

dz 

dt 


[t] 

[0] 

[0] 


<le_  ( 

Erf  — 

me  V 

Qe_  ( 

—Edc 

me  V 

Vq  cos 

[0] 

»-4) 

o  dy 
dt 


>RF 


—Vo  sin  [0 


When  solved  these  equations  yield  the  terms 


(30) 

(31) 

(32) 

(33) 


vz[t] 


z[i] 


Edc  ~  Edc  cos 


BRFQet 

rrte 


+  Brfv  0  cos 


BRFQet 


me 


+ 


Erf  +  sin 


^RFlet 


Erf  —  Erf  cos 


BRFqet 

me 


+  Edc  sin 


Brf 

BRFqet 

me 


Brfv  0  sin 


BRF<let 


+ 


(35) 


— Brf  (  EDCqet  -  2 me  sin 

BRFQe  V 


Brf 
BRFqet 
2  mP 


Edc  cos 


BRFqet 
2  mP 


-Brfv o  cos 
1 


BRFqet 
2  me 


+ 


BRFQe  V 
+BRFmev0  cos 


(  EDCme  -  EDCme  cos 
BRFqet 


+  Erf  sin 
BRFqet 


rnf 


BRFqet 
2  me 

-  BRFmev o  cos  I 


rrip 


+ 


+  BRFERFqet  -  ERFme  sin 


BRFqet 


nip 


(36) 


(37) 


The  value  for  z[t\  may  be  numerically  solved  to  get  the  particle  transit  time,  which 
using  the  same  process  used  to  generate  Equation  24  allows  for  the  calculation  of  the 
energy  associated  with  the  particle  at  its  time  of  impact  with  the  dielectric  window. 

However,  a  further  simplification  of  these  equations  is  obtained  when  the  magnetic 
field  influence  is  negligible.  In  the  limit  of  BRF  approaches  zero,  the  equations  of 
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motion  indicated  above  simplify  to  the  form  of 


Vy[t] 

=  v0  cos  [<f>]  +  —ERFt 
rne 

(38) 

Vg[t] 

=  -v0  sin  [(f)]  d  —EDCt 
rne 

(39) 

y[t] 

=  Vo  cos  [4>\  t  d  \ — ERFt2 

2  me 

(40) 

z[t ] 

=  -v0  cos  [4>]t  d  \—EDCt 2 

2  me 

(41) 

Where  the  transit  time  r  is  Equation  15.  For  a  normal  ejection  (0  =  |)  the  final 
energy  gained  during  the  particle’s  flight  is 


£ Strike  m  d 


me(\E2RFvl  +  \E2DCvl  +  | ErfVq) 

(Edc 


(42) 


After  replacing  the  final  calculation  of  the  boundary  curve  equations  and  associating 
initial  velocity  with  initial  energy  the  ERF  held  may  be  written  as 

4  \Z^0m  Strike  ^Om)  2(£ Strike  ^Om)  Sm  [r/l] 

£ Strike  d  (3^0m  d  £ Strike )  COS  (0] 

showing  the  linear  relationship  that  one  would  anticipate  for  the  boundary  curves  in 
the  static  held  case. 


2.2  Method 

As  previously  stated,  when  space  charge  effects  are  omitted  modeled  particles  have 
no  associated  fields  and  when  emitted  from  the  dielectric  do  not  have  leave  behind  a 
positive  charge  and  a  bias  held  must  be  applied  in  the  simulation  in  order  to  attract 
particles  back  to  the  surface  of  the  dielectric  enabling  multipactor  growth.  In  the 
NSC  case,  the  simulation  results  may  also  be  compared  to  analytic  solutions  for  the 
boundary  curves  allowing  for  additional  verification  of  computations  run  in  ICEPIC. 
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Already  both  Ficthl  [8]  and  Rogers  [19]  have  had  great  successes  in  modeling  lower 
boundary  curves  comparable  to  those  of  Kishek  [14]  and  Lau  [13]. 

Defining  the  EM  Radiation. 

Within  the  ICEPIC  code,  there  is  a  capability  to  create  a  plane  wave  shape  within 
the  simulation  space.  This  plane  wave  function  allows  for  the  simulation  of  electro¬ 
magnetic  radiation  expected  at  the  waveguide  center.  The  plane  wave  modeled  the 
electromagnetic  radiation  associated  with  a  TEW  mode  propagating  approximately 
at  the  waveguide  center.  The  electric  field  amplitude  of  the  plane  wave  could  be 
defined  by  the  parameter  EO  and  the  frequency  could  be  defined  with  the  frequency 
parameter.  A  picture  of  the  plane  wave  propagation  through  the  dielectric  may  be 
seen  in  Figure  8.  Note  that  PML’s  were  used  to  prevent  backscatter  of  the  plane 
wave  from  the  edges  of  the  simulation.  The  plane  wave  shape  defines  the  extent  of 
the  propagating  electromagnetic  radiation,  a  second  dielectric  was  also  used  to  pre¬ 
vent  the  back  scatter  of  particles  from  interacting  with  the  edges  of  the  plane  wave 
shape. 

Particle  Emission. 

In  order  to  define  secondary  particle  emission  special  particle  emissions  must  be 
defined  to  encompass  any  shape  from  which  particle  emission  was  desired.  In  order 
to  understand  the  process  of  particle  emission  in  ICEPIC,  it  becomes  important  to 
take  a  snap  shot  of  the  ICEPIC  manual’s  definition  of  particle  emission  [5].  Since 
Ficthl  and  Rogers  simulations  [8]  [19],  the  interactions  between  emitted  BEAM  and 
SECONDARY  particle  emission  have  been  updated  to  classify  secondary  emission.  The 
parameter  second_species_name  defines  the  species  with  which  a  particular  particle 
region  interacts.  This  improvement  results  in  three  different  particle  shapes  being 


22 


Figure  8.  The  plane  wave  shape  defined  within  the  ICEPIC  environment  along  with 
the  expected  region  of  multipactor  within  the  system  setup. 

used  to  define  primary  emitted  particles,  the  secondary  emitting  region  that  reacts 
to  the  presence  of  primary  particles,  and  a  secondary  emitting  region  that  reacts  to 
secondary  particles  emitted.  In  the  ICEPIC  code,  the  primary  emitted  particles  are 
defined  by  direction  injection,  the  rate  of  injection  and  their  average  temperatures. 

Secondary  emission  shapes  have  a  few  different  properties  than  that  associated 
with  primary  emission  regions.  They  follow  the  empirical  equations  of  secondary  elec¬ 
tron  emission  developed  by  Vaughan  [24] .  Material  parameters  out  of  these  equations 
are  used  to  define  secondary  emission  from  this  particle  shape,  as  well  as  scattering 
and  reflection  coefficients  dealing  with  the  material’s  interaction  with  the  impacting 
particles.  The  parameters  used  for  the  mapping  of  the  boundary  curves  followed  from 
values  used  by  Fictlil,  Rogers,  and  Kishek  [8]  [19]  [15]  and  are  reproduced  in  Table  1. 
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Table  1.  Parameters  from  Ficthl  [8]  used  to  verify  ICEPIC  simulation  setup,  and  in 
remapping  the  susceptibility  diagrams. 


Umax  0 

2 

C 

'-'max  0 

400  eV 

ks 

1 

T 

2.1  eV 

Resolution 

3000 

Static  Field  Implementation. 

For  the  NSC  case,  a  bias  field  is  required  to  return  particles  to  the  surface  within 
the  simulation.  Within  the  simulation  environment  static  uniform  electric  fields  are 
able  to  be  defined  for  the  entire  simulation  environment  with  the  default  _e  command. 
Uniform  magnetic  fields  may  also  be  assigned  to  the  environment  with  the  def  ault_b 
command.  By  applying  static  electric  field  within  the  simulation  in  the  —z  direction, 
the  bias  field  associated  with  the  NSC  mapping  may  be  achieved. 

Building  Susceptibility  Curves. 

By  parametrically  varying  the  values  for  the  plane  wave  field  amplitude  and  the 
bias  electric  field,  the  susceptibility  diagram  was  constructed.  As  seen  in  Figure  9  and 
Figure  10,  the  construction  of  the  boundary  curves  involved  mapping  the  outcome, 
growth  or  decay,  resulting  from  a  survey  of  E\\  values  for  a  single  fixed  value  of  the 
bias  electric  field,  E± .  This  was  then  repeated  for  a  sequence  of  new  bias  electric 
field  inputs.  The  resulting  points  differentiating  growth  and  decay  regions  were 
then  fitted  with  a  linear  regression  corresponding  to  the  linear  relationship  observed 
in  Equation  25.  The  slopes  were  then  recorded  and  compared  to  theory. 
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E\  ( MV|  m) 

Figure  9.  Boundary  curve  generated  from  observing  the  growth  and  decay  of  multi- 
pactor  by  parametrically  surveying  E±  and  E\\  for  given  parameters.  The  numbers  in 
bubbles  correspond  to  growth  and  decay  particle  populations  represented  in  Figure  10. 

2.3  Results  and  Analysis 

Previous  results  in  characterizing  boundary  curves  in  the  susceptibility  diagram 
were  reported  by  Kishek  and  Lau  [15].  These  follow  the  linear  boundary  conditions 
indicated  in  Figure  11,  which  confirmed  the  results  previously  indicated  for  the  theo¬ 
retical  boundary  curve  indicated  in  Equation  25.  In  Fichtl’s  and  Roger’s  work  only 
the  lower  boundary  curve  was  analyzed  to  confirm  the  validity  of  ICEPIC  measure¬ 
ments  [8]  [19].  However  when  the  boundary  curve  associated  with  the  high  crossover 
energy  (e2 )  was  analyzed  in  ICEPIC,  a  large  discrepancy  was  observed  between  the 
remapping  and  previous  results  from  Monte  Carlo  and  other  PIC  code  simulations 
[3]  [13]  [15]  [22]  [12],  The  new  susceptibility  mapping  is  shown  in  Figure  12  .  After 
examining  the  simulation,  the  source  of  this  discrepancy  was  traced  to  the  use  of  the 
plane  wave  feature  in  ICEPIC.  In  previous  techniques  for  mapping  susceptibility  di- 


25 


(a)  Well  above  boundary  curve  at  with  field  val¬ 
ues  of  E\\  =  3.3 MV/m  and  E±  =  1  MV/m 


Time  (s)  *  10  » 


(c)  Just  above  boundary  curve  at  with  field  val¬ 
ues  of  E\\  =  2.5 MV/m  and  E±  =  1  MV/m 


(b)  Well  below  boundary  curve  at  with  field  val¬ 
ues  of  E\\  =  1.8 MV/m  and  E±  =  1  MV/m 


(d)  Just  below  boundary  curve  at  with  field  val¬ 
ues  of  E\\  =  2.3 MV/m  and  E±  =  1  MV/m 


Figure  10.  A  temporal  variation  of  total  particle  population  in  the  vicinity  of  growth 
and  decay  around  the  boundary  curve.  The  numbers  in  the  bubbles  correspond  to  the 
numbers  in  Figure  9  showing  the  process  to  create  a  susceptibility  diagram. 

agrams,  Monte  Carlo  techniques  have  only  considered  the  electric  field  component  of 
electromagnetic  radiation  propagating  through  the  waveguide  dielectric.  In  contrast, 
the  plane  wave  shape  in  ICEPIC  also  contains  the  magnetic  field  component  present 
in  the  radiation  propagating  through  the  dielectric. 

In  order  to  confirm  the  magnetic  field  as  the  origin  of  the  discrepancy,  another 
mapping  was  accomplished  using  static  electric  and  magnetic  fields  to  represent  the 
incoming  EM  radiation.  These  results  did  confirm  the  origin  of  the  shift  of  the  curve 
associated  with  higher  energy,  e2 ,  boundary.  This  shift  is  shown  in  Figure  13,  where 
the  change  in  the  slope  of  the  boundary  curve  associate  with  e2is  apparent. 

To  summarize,  the  lower  boundary  curve  associated  with  the  first  crossover  energy 
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£±(MV/m) 


Figure  11.  Recreation  of  previous  work  by  Kishek  and  Lau  [15].  The  boundary  curves 
have  been  fitted  with  a  linear  regression  to  capture  their  extent. 

£i ,  tends  to  follow  closely  with  previous  work  done  on  analyzing  the  boundary  curves. 
The  slope  of  this  first  crossover  energy  is  slightly  higher  in  the  Estatic  case.  However, 
in  the  case  of  the  upper  boundary  curve  associated  with  the  impact  energy  e2  the 
slope  comparison  confirms  that  the  magnetic  held  causes  multipactor  growth  at  high 
enough  E\\  values  despite  the  presence  of  no  return  held,  E±  =  0. 

These  differences  in  slope  change  between  the  upper  and  lower  boundary  curves 
can  clarified  by  considering  the  Lorentz  force  imparted  by  the  magnetic  held.  By 
again  using  the  linear  slope  analogy,  Equation  25  and  Equation  43  lead  to  the  form 
of  the  boundary  curve 

E\\  =  7-Ej_eif  (44) 

where  7  is  defined  as  the  slope  of  the  boundary  curve  produced  with  consideration 
of  only  the  electric  held  and  E±eff  is  the  effective  electric  held.  By  defining  E±eff  as 
a  function  of  the  magnetic  and  electric  held  for  a  particle  with  an  average  velocity 
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Figure  12.  Susceptibility  diagram  results  based  on  the  ICEPIC  outputs  using  the  plane 
wave  shape  command.  Includes  both  a  dynamically  varying  Ey  and  fields. 


<  vavg  >  for  the  length  of  transit  leads  to  the  equation 


E\\  —  7 (EDC+  <  V  >avg  B |j ) 


and  considering  B\ 


obtain 


E\ 

E\ 


7  (Edc  + 

7  Edc 
(1-7^ 


<  v  > 


avg 


E,\) 


which  we  may  relabel  as  a  new  effective  slope  field  to  form 


(45) 


(46) 


where 


E\\  —  7 BfieldEoC 


(47) 


7 


(1-7 


<V>a 


'JB  field 


Figure  13.  The  results  of  mapping  the  boundary  curves  for  the  static  field  case.  Note 
the  shift  that  occurs  when  a  magnetic  field  whose  amplitude  is  equivalent  to  the  am¬ 
plitude  of  a  magnetic  field  in  a  plane  wave  (13  =  ^-). 

This  equation  establishes  the  shift  in  the  curves  based  on  the  inclusion  of  the  magnetic 
held  and  leads  to  the  results  pictured  in  Figure  14.  As  one  can  see,  there  is  very  close 
agreement  to  the  expected  shift  actually  observed  in  the  remapping  of  the  magnetic 
held  into  a  new  form.  A  summary  of  linear  regression  results  are  presented  in  Table  2. 

The  remapped  susceptibility  diagram  indicates  an  opportunity  to  mitigate  the 
initiation  of  multipactor  with  the  application  of  an  external  bias  held.  From  Figure  12, 


Table  2.  Summary  of  linear  regressions  for  boundary  curves  illustrated  in  figures  11, 
12,  13.  The  values  are  organized  by  their  association  with  the  lower  (ei  )  and  higher 
(e2  )  cutoff  energies  associated  with  electron  impact. 


Study 

Boundary  Curve  £1 

Boundary  Curve  e2 

Slope 

Intercept 

Slope 

Intercept 

Original  Kishek  [15] 

1.54 

0.1 

28.1 

1.2 

Planewave 

2.01 

0.43 

-34.2 

0.48 

E static 

2.22 

-0.03 

57.67 

0.10 

E static  ckQd.  B static 

2.25 

-0.04 

-20.39 

-.53 
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Lower  Energy  Boundary  Curve 

Higher  Energy  Boundary  Curve  Including  Magnetic  Field 
Higher  Energy  Boundary  Curve  Excluding  Magnetic  Field  I 
Theoretically  Predicted  Upper  Boundary  Curve 


£i(MV/m) 


Figure  14.  Theoretical  line  fitting  of  the  shifted  upper  boundary  in  the  static  field 
based  upon  Equation  47  showing  agreement  between  in  the  shift  observed  with  the 
theoretical  model  including  the  magnetic  Lorentz  force. 


the  required  bias  to  prevent  secondary  emission  growth  is  much  less  if  it  is  applied 
in  the  —  z  direction  in  respect  to  the  simulation  setup.  This  has  been  previously 
indicated  in  Spaulding’s  work  in  characterizing  the  application  of  an  external  bias  field 
in  preventing  the  growth  of  nmltipactor  [21]  and  is  confirmed  in  the  new  susceptibility 
mapping. 
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III.  Analysis  and  Simulation  of  Reflection 


The  multipactor  mitigation  problem  stems  from  the  conventional  properties  of  a 
plasma’s  interaction  with  high  intensity  electromagnetic  energy  at  microwave  frequen¬ 
cies.  These  frequencies  typically  range  from  1-3  GHz  and  tend  to  lead  the  electron 
cascade  process  at  1-3  MV/m  electric  field  intensities.  At  these  levels  the  plasma 
will  begin  to  exhibit  some  interesting  effects,  primarily  the  reflection  of  microwave 
radiation  at  the  dielectric  window  and  the  heating  of  the  dielectric  window.  In  this 
chapter,  the  reflection  associated  with  multipactor  plasma  is  characterized  with  both 
analytic  development  and  simulation.  In  order  to  properly  characterize  the  electron 
cloud  formed  along  the  dielectric,  the  Space  Charge  must  be  included  in  the  particle’s 
interaction  with  the  electromagnetic  field  and  with  each  other. 

3.1  Theory 

Understanding  Multipactor  in  Real  Time. 

In  the  SC  case,  the  incoming  microwave  radiation  is  in  the  form  of  a  TE\$  mode 
propagating  inside  the  waveguide  before  propagating  through  the  dielectric  surface. 
Previous  experiments  completed  by  Neuber  at  Texas  Tech  [16]  indicate  that  the  region 
most  prone  to  multipactor  is  along  the  center  of  the  waveguide  window  corresponding 
to  the  greatest  values  for  E\\  in  the  TUlj0  mode.  In  this  region  of  the  electromagnetic 
mode,  the  microwave  radiation  may  be  approximated  with  as  a  plane  wave  normally 
incident  on  the  dielectric  window. 

An  example  of  particle  growth  associated  with  the  space  charge  case  is  shown  in 
Figure  15.  It  illustrates  that  the  particle  growth  follows  a  linearly  increasing  oscilla¬ 
tion,  with  the  frequency  of  oscillation  in  particle  numbers  corresponding  to  twice  the 
frequency  of  the  propagating  electromagnetic  radiation.  In  the  SC  case,  particles  are 
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Figure  15.  The  dynamic  particle  growth  associated  in  a  multipacting  system.  The 
propagating  electromagnetic  energy  causing  the  particle  growth  has  a  held  amplitude 
of  E\\  =  3MV/m  and  a  frequency  of  /  =  2A5GHz.  Every  two  periods  of  particle  oscillation 
corresponds  to  a  single  period  of  the  microwave  radiation. 


returned  to  the  surface  of  the  dielectric  primarily  by  compensating  positive  surface 
charge  left  on  emitting  surface.  By  parametrically  plotting  the  E\\  held  associated 
with  the  plane  wave  and  the  E±  held  associated  with  the  surface  charge  oscillation,  a 
“bowtie”  curve  is  constructed.  The  “bowtie”  curve  is  a  valuable  tool  in  characterizing 
the  temporal  variation  of  multipactor. 

By  following  the  arrows  in  the  “bowtie”  curve  plotted  in  Figure  16,  the  evolution 
of  secondary  particle  emission  begins  with  the  E\\  held  supplying  energy  to  shift  the 
average  impact  energy  of  electrons  into  a  region  of  multipactor  initiation.  These  par¬ 
ticles  continue  to  impact  with  the  surface  causing  an  electron  avalanche  that  leads  to 
an  increasing  E±  held.  Eventually,  the  E\\  held  begins  to  decrease  corresponding  with 
the  sinusoidal  oscillation  of  the  plane  wave  propagating  through  the  dielectric.  This 
leads  to  lowering  the  average  impact  energy  of  particles  in  the  electron  cloud.  Eventu¬ 
ally,  the  average  impact  energy  falls  below  the  boundary  curve,  leading  to  multipactor 
decay  and  a  consequent  decrease  in  the  number  of  particles  present  in  the  electron 
cloud.  This  process  is  than  repeated  for  negative  values  of  the  E»  held,  indicating  that 
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Figure  16.  “Bowtie”  curves  based  on  Kim  and  Verboncoeur’s  previous  work  [12]  in 
characterizing  multipactor  growth  and  decay.  This  curve  represents  the  steady  state 
dynamics  of  the  E\\  and  E±  fields  parametrically  plotted  in  relation  to  each  other.  The 
red  lines  correspond  to  boundary  curves  associated  with  separating  regions  of  multi¬ 
pactor  growth  and  decay.  The  arrows  shows  the  direction  that  the  fields  evolve 


only  the  field  amplitude  is  important  in  determining  multipactor  growth  and  decay. 
It  is  important  to  understand  this  evolution  of  particle  count  when  characterizing  the 
reflection  of  electromagnetic  radiation. 


Calculating  the  Power  Transmitted  through  the  Multipactor. 

The  reflection  problem  of  microwave  radiation  stems  from  the  underlying  issues 
concerning  all  plasmas  reaction  with  electromagnetic  radiation.  As  seen  in  Fox’s 
textbook  [9],  the  theory  behind  bulk  plasmon  effects  in  materials  could  be  applied  to 
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the  microwave  multipactor  problem.  Starting  with  Maxwell’s  equations 


V  •  E  = 

P_ 

(48) 

eo 

VB  = 

0 

(49) 

dB 

(50) 

V  x  E  = 

~~dt 

d  E 

1 

(51) 

j  +  £°~dt 

—V  x  B. 
ho 

one  may  develop  the  wave  equation  for  electromagnetic  radiation  by  taking  the  deriva¬ 
tive  of  equation  51  and  substituting  Equation  50  [9] 


dj  d2E 
dt+e°^E 


V  x  (V  x  E). 
ho 


(52) 


The  electrons  will  move  according  to  the  electric  field  mev  =  —qeE  and  current  may 
be  expressed  as  j  =  —neqev.  Using  these  two  equations  allows  one  to  write  the  change 
in  current  density  as  %  =  ^^-E.  Substituting  this  and  c2  =  — —  into  Equation  52 
reveals: 

^  +  (M i)E  =  — c2V  x(VxE)  (53) 

otz  eo  me 

2 

where  is  normally  identified  the  plasma  frequency: 

eo  me  J  r  1  J 


UJ2P  = 


Neq , 


e0me 


(54) 


Knowing  that  the  electric  field  may  be  split  into  its  transverse  and  longitudinal  com¬ 
ponents  allows  one  to  write  equation  53  as  two  independent  field  equations: 


d2Et 

dt2 


+  u2Et  -  c2V2Et  =  0 


(55) 
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and 


^+^4  =  0.  (56) 

Assuming  a  plane  wave  solution  for  Equation  55  the  electric  field  of  the  traveling 
wave  is 

E\\  =  E0e (57) 
where  the  relation  between  k  and  u  is 


(58) 


Observe  that  radiation  with  an  angular  frequency  less  than  the  plasma  frequency 
will  not  propagate  through  the  extent  of  the  plasma  [9].  With  minor  rewriting  of 
Equation  58  the  refractive  index  for  the  plasma  may  then  be  found. 


n 


2 


1  - 


(59) 


Assuming  propagation  in  the  ^-direction  and  writing  k  in  terms  of  the  refractive  index 
n,  we  obtain  a  new  form  of  Equation  57 


E\ 


(60) 


Notice  that  if  the  refractive  index  and  the  wave  vector  become  pure  imaginary,  the 
field  amplitude  will  decay  exponentially.  This  will  occur  if 


Up  >  u. 


(61) 
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The  transverse  field  solution  shown  in  Equation  60  may  then  be  rewritten  as 


B,  =  E0e~f“ M* 


(62) 


where 

(\A>p(2)2  “  ^2) 

a(z )  =  ^ 

c 

By  using  Valfell’s  [22]  particle  density  distribution  characterizing  the  extent  of  the 
multipactor  cloud, 


where 


and 


n(z ) 


n0 

(1  +  ^)2 


(64) 


(65) 


Equation  63  may  be  rewritten  as 


a(z)  = 


This  can  be  further  simplified  by  writing 


where 


a(z) 


i  (yi-z^w2) 

A0  p(z) 


P(z)  = 


CO  CO  z 

- 1 - • 

UJp  vt 


(66) 


(67) 
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Integrating  f  a(z)  dz  over  the  range  of  the  multipactor  extent  leads  to 
J a(z)  dz  = 


Equation  68  leads  to  a  wave  equation  for  the  electric  held  transmitted  through  the 
full  extent  of  the  multipactor  plasma 


a(P)~rkd/3  —  —  f  a(/3)d/3 
dp  u  J 


vt 


(v/rqp) 


u  £ 


P 


d/3 


Ln[£]  +  Ln  [l  +  y/l  -  ? 


(68) 


t/j  =  E0e“?(_^-LnK1+Ln[i+^])e-u 
=  /.'i£  '  f  V  1.'  [l+V  i-£2])  e-iwt 


(69) 


Knowing  that  the  intensity  of  the  microwave  radiation  related  to  the  electric  held  is 
a  function  of  the  Poynting  vector  [10] 

\S\  =  e0c^2,  (70) 

allows  for  Equation  69  to  be  rewritten  as  a  function  of  intensity 

|5Tr„|  =  £0e£0V?(^+M«-.4+^]).  (71) 
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Equation  71  may  be  then  converted  to  the  transmitted  power  by  multiplying  the 
intensity  by  the  area  of  the  plane  wave 


where 


R 


Transmitted 


=  ( AxAy)e0cE( 


1  +  A/7 1 


£  =  -■ 
(jjqr) 


(72) 


Rewriting  Equation  72  as  a  new  function  of  — 


where 


(73) 


Figure  17  illustrates  the  curve  relating  the  transmission  coefficient  to  the  ratio  of 
the  plasma  frequency  to  the  angular  frequency  of  the  propagating  electromagnetic 
radiation.  The  transmission  coefficient  is  based  on  the  angular  frequency  of  the 
propagating  EM  radiation  (a;),  the  emitted  electron  velocity  that  is  dependent  upon 
material  constituting  the  dielectric  (vt),  and  the  plasma  frequency  based  upon  the 
particle  density  at  the  dielectric  surface  (c op).  The  only  value  reliant  on  the  electric 
held  amplitude  of  the  propagating  microwave  radiation  is  going  to  be  the  plasma 
frequency  due  to  it’s  dependence  on  the  particle  density.  As  particle  count  increases 
the  transmission  of  radiation  decreases  through  the  multipactor.  If  the  particle  den¬ 
sity  becomes  great  enough,  there  will  be  no  transmission  through  the  multipactor 
approaching  the  limiting  case  of  a  total  plasma  reflection  associated  with  plasmas 
where  the  plasma  extent  is  much  larger  than  the  wavelength  of  the  propagating  elec- 
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Figure  17.  The  curve  relating  transmission  coefficient  to  the  ratio  of  the  plasma  fre¬ 
quency  (uip)  to  the  angular  frequency  oj  of  the  plane  wave  propagating  through  the 
multipactor  generated  from  Equation  73.  The  electron  velocity  ( vt )  corresponds  to  an 
electron  emitted  with  an  energy  of  =  2eV. 


tromagnetic  radiation. 


Calculating  the  Intensity  Reflected. 

In  order  to  calculate  the  power  per  unit  area  reflected  for  microwave  radiation 
incident  on  the  dielectric  surface,  one  may  relate  Equation  72  to  the  incident  plane 
wave  intensity  by 


PRe  fleeted  Plncident  (PDeposited  “t"  pTransmitted'}  ■ 


(74) 


In  which  PoeposUed  is  defined  by  Valfells  [22]  to  be 
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with  this  final  equation  the  reflected  power  may  be  calculated  analytically  from 
Valfell’s  equation  for  power  deposition  and  the  new  equation  developed  for  the  inten¬ 
sity  of  the  electromagnetic  radiation  transmitted  through  the  multipactor. 

In  order  to  relate  the  transmitted  power  in  Equation  72  and  the  power  deposited 
from  Equation  75  to  the  number  of  particles  indicated  in  the  ICEPIC  output  files,  the 
charge  density  at  the  surface  of  the  dielectric,  n0,  had  to  be  related  to  the  total  charge 
contained  in  the  multipactor  cloud  iVrot •  This  may  be  accomplished  by  expressing 
iVTot  as  a  function  of  n0  using  Valfcll’s  charge  distribution  [22] 


iVrot  = 


=  vt 


poo  p  Ay  p  Ax 

/  /  7 
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l  i  <^p±.  \  2 
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mee0 


A  A  /  n0me^0  *  a 

AxAy  =  vH  x — AxAy. 


de 


(77) 


With  this  derived  relationship  the  total  power  transmitted  through  the  multipactor 
may  be  calculated  from  the  total  amount  of  particles  in  the  multipactor  (IVrot)  • 


3.2  Methodology 

Reflection  of  the  Dielectric  Window. 

Before  quantifying  the  reflection  of  multipactor  within  the  cavity,  the  reflection 
associated  with  the  dielectric  window  must  be  considered  in  order  to  properly  identify 
reflection  associated  with  the  dielectric  and  reflection  associated  with  the  multipactor. 
Within  ICEPIC,  the  dielectric  constant  of  the  windows  is  set  to  be  ewindow  to  1.0001, 
resembling  the  dielectric  constant  associated  with  anti-reflection  coatings  present  on 
current  window  dielectrics  used  in  HPM  devices  [4],  The  actual  reflection  coefficient 
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is  calculated  to  be 


R  =  |  — |2  «  |^^|2  =  6.24  x  1(T10,  (78) 

ti  +  1  y  e  +  1 

indicating  that  the  reflection  at  the  dielectric  windows  for  the  simulated  value  of  e 
is  negligible  in  comparison  to  the  reflection  value  of  the  multipactor  found  by  Ficthl 
[8],  justifying  it  being  ignored  within  this  analysis. 

Using  PROBE  to  Determine  Electric  Fields. 

The  determination  of  electric  fields  associated  with  reflection  and  transmission 
was  completed  using  the  PROBE  dump  command  in  ICEPIC.  This  was  a  useful  tool 
in  analyzing  the  electric  field  at  specific  points  within  the  simulation.  Each  selected 
z  location  was  given  seven  (matching  Nx)  uniformly  spaced  PROBE  dumps  creating 
a  probe  line  which  when  averaged  allows  for  the  elimination  of  noise  in  the  analysis 
of  the  electric  and  magnetic  fields.  Locations  of  these  probe  lines  are  illustrated  in 
Figure  18.  Probe  lines  were  chosen  primarily  to  analyze  two  different  parts  of  the 
system,  the  reflection  of  the  planar  wave  and  the  parallel  and  perpendicular  electric 
fields  associated  with  the  system,  E^ef  measures  the  electric  field  amplitude  reflected 
from  the  particle  cloud  in  front  of  the  dielectric  window,  ETrans  corresponds  to  the 
field  transmitted  through  the  particle  cloud,  Einc  is  used  to  measure  the  field  initially 
launched  by  the  [Planewave]  section  of  the  code,  and  Edc  is  used  to  measure  fields 
perpendicular  and  parallel  to  the  surface  of  the  dielectric.  However,  the  Ejnc  probe 
will  contain  a  superposition  of  the  reflected  and  incident  waves,  as  shown  in  equation 

Einc  =  Eorig  +  E^ef  (79) 
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Figure  18.  Various  “probe  line”  locations  within  the  simulation.  They  indicate  the 
position  of  the  probes  measuring  Eftef  to  be  in  a  place  that  will  receive  only  the  reflected 
EM  radiation  generated  by  the  plane  wave.  The  probes  measuring  Ejnc  measure  the 
initial  generated  electric  field.  Probe  Ejjc  measures  the  electric  fields  close  to  the 
dielectric  to  ascertain  the  relationship  between  E±  and  E\\  .  The  probes  known  as  Exrans 
measure  the  EM  radiation  transmitted  through  the  dielectric. 


where  E0rig  is  the  original  emitted  electric  held.  As  indicated  by  equation  78  there 
will  be  very  little  reflection  from  the  dielectric  windows  in  the  model.  Therefore, 
within  the  simulation  reflected  electromagnetic  energy  will  be  primarily  attributed 
to  the  multipactor  on  the  outside  of  the  dielectric  window.  The  rest  of  the  material 
properties,  were  based  on  Verboncoeurs  [12]  and  Rogers  [19]  previous  work  for  char¬ 
acterizing  secondary  electron  emission  contained  in  Table  3.  With  this  knowledge  of 
held  oscillations  received  from  the  PROBE  dump,  a  simple  calculation  [10]  using 


{S)  = 

f  <£>2 

(80) 

< p )  = 

(S')XSySySyS 

(81) 
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Table  3.  Parameters  from  Rogers  [19]  used  to  verify  ICEPIC  simulation  setup. 


OmaxO 

3 

£max  0 

420  eV 

ks 

1 

T 

2.1  eV 

Resolution 

12000 

(a)  Method 


Figure  19.  Using  the  results  from  the  Probe  Ejjc  to  manufacture  the  “Bowtie”  curves 
associated  with  varying  electric  fields  for  an  electric  field  amplitude  of  E\\  =  3MV/m 
and  /  =  2.45GHz.  " 

where  xsys  and  ysys  are  the  respective  dimensions  of  the  simulation  in  the  x  and  y 
directions  enables  the  calculation  of  power  being  transmitted  as  a  function  of  the 
electric  held  in  the  system.  As  indicated  in  Section  3.1  the  “bowtie”  curves  are 
excellent  tools  in  determining  the  characteristics  associated  with  the  temporal  vari¬ 
ation  of  the  multipactor  plasma.  By  measuring  both  the  E±  and  E»  fields  with  the 
probe  Eoc  line,  the  parametric  relationship  between  both  En  and  E±  may  be  plotted 
as  shown  in  Figure  19. 
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(a)  Method 


(b)  Typical  Output 


Figure  20.  Using  Equation  82  to  derive  the  E±  field  based  on  the  particle  count  and 
weight  present  in  the  multipactor  cloud. 

Using  Particle  Count  to  Determine  Electric  Field  Perpendicular. 

The  benefit  of  knowing  the  particle  count  in  the  space  charge  case  comes  from  the 
ability  to  calculate  the  E±  held  from  the  charged  particles.  If  we  treat  the  sheet  of 
positive  charge  on  the  dielectric  as  being  a  surface  charge,  and  also  treat  the  electric 
held  produced  by  the  electrons  in  the  same  manner,  the  superposition  of  the  electric 
helds  leads  to  the  equation 

I  •  Q  T  ot 

6qA 

where  Qrot  is  defined  as  the  total  charge  associated  with  the  window  dielectric  and 
A  is  the  area  of  the  simulation  space  in  the  x-y  plane.  A  “bowtie”  trajectory 
constructed  of  using  the  particle  count  to  generate  the  E±  held  is  displayed  in  Fig¬ 
ure  20. This  technique  also  removes  the  noise  associated  with  the  method  of  using  the 
probe  Edc  line  measurements  pictured  in  Figure  19.  The  methods  show  that  there  is 
a  discrepancy  that  exists  between  the  values  reported  by  both  methods. 
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Particle  Weighting. 


Due  to  the  requirement  for  computational  efficiency,  the  ICEPIC  code  does  not 
treat  electrons  as  single  particles  in  the  system.  Rather,  numbers  of  electrons  are 
grouped  together  into  single  pseudoparticles.  It  is  therefore  necessary  to  establish  the 
total  charge  available  in  the  multipactor  cloud  within  an  accurate  particle  weighting 
to  relate  the  number  of  macroparticles  in  the  ICEPIC  code  to  the  number  of  electrons 
being  represented  by  each  macroparticle.  This  becomes  important  for  the  calculation 
of  both  the  reflection  and  the  evaluation  of  the  E±  field. 

The  particle  weighting  is  initially  determined  in  the  simulation  based  on  the  beam 
parameters  associated  with  the  initial  current  injected  in  the  simulation  to  initiate 
multipactor.  The  calculation  of  this  weighting  is  based  upon  the  the  charge  associated 
with  the  macroparticle  [5] 


Qinacro 

which  combined  with 
weight 


current  inject  .interval  dt 
inject 


A  Cell 
Emission 


Q.macro 


Qe 


(83) 


(84) 


Definition  of  these  quantities  used  in  Equation  83  are 

•  current  -  This  dictates  an  overall  current  (  in  Amps)  maintained  by  the 
emitting  surface.  This  plays  into  the  weighting  the  weighting  of  particles  by 
dictating  a  required  charge  per  shape  injected  into  the  system. 

•  inj  ect  -  This  is  the  amount  of  particles  injected  from  each  cell  of  the  indicated 
injection  shapes.  As  one  would  suspect  this  also  plays  into  the  weighting  of 
particles  in  the  multipactor  simulation. 

•  inj  ect  .interval  -  This  is  the  time  in  terms  of  dt  between  ejections  from 
each  cell.  An  important  feature  to  note  is  if  there  are  less  particles  than  cells 
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Table  4.  Comparison  of  weighting  in  terms  as  outputted  by  the  ‘ice.stat’  file,  the  PART 
dump  file  and  the  theoretical  output  of  equation  84  for  the  two  dimensional  NSC  case 
of  simulations,  where  the  baseline  value  of  the  the  cases  has  a  current  of  1,  a  Nx  of  3, 
an  inject  of  1,  an  inject  .interval  of  1  and  a  resolution  of  12000. 


Property  Changed 

New  Value 

Particle  Dump  (C) 

‘ice.stat’  (C) 

Theoretical  (C) 

Baseline 

7.94e-15 

7.94e-17 

7.94e-15 

current 

100 

7.94e-13 

7.94e-15 

7.94e-13 

Nx 

50 

4.76e-16 

4.76e-18 

4.76e-16 

inject 

7 

1.13e-15 

1.13e-17 

1.13e-15 

inj  ect.interval 

5 

3.97e-14 

3.97e-16 

3.97e-14 

Resolution 

3000 

3.17e-14 

3.17e-16 

3.17e-14 

being  injected  at  a  time  interval,  the  cell  faces  will  randomly  eject  particles. 

•  dt  -  The  timestep  of  the  simulation  defined  by  Equation  1  in  Chapter  I 

•  Aceii  -  The  area  associated  with  a  simulation  unit  grid  (defined  as  dz 2). 

•  A  Emission  -  The  area  of  the  emitting  particle  phase,  defined  as  the  dimension 
in  the  x  and  y  directions. 

Effective  use  of  the  Equations  83  and  84  enables  efficient  usage  of  computational 
resources  within  the  simulation  environment.  However,  care  must  be  taken  to  balance 
the  need  for  computational  efficiency  and  the  number  of  particles  required  to  maintain 
statistical  validity  of  the  ICEPIC  simulation  . 

When  verifying  Equations  83  and  84  a  discrepancy  was  noticed  in  the  values  of 
the  charge  associated  with  the  macroparticle  in  the  ‘ice.stat’  hie  utilized  for  tracking 
particle  count.  By  varying  the  parameters  inside  Equation  83  it  was  determined  that 
the  ‘ice.stat’  hie  always  outputted  macroparticle  charges  a  factor  of  a  hundred  less 
than  those  predicted.  These  inconsistencies  reported  for  the  NSC  case  may  be  viewed 
in  Table  4.  Further,  investigation  confirmed  that  this  inconsistency  also  occurs  within 
the  SC  case.  Both  of  these  cases  were  ran  for  a  two  dimensional  simulation.  The  PART 
dump  hies  that  describe  particle  data  were  also  investigated  for  the  particle  charge 
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reported  and  as  shown  in  Table  4,  concurred  with  Equations  83  and  84.  Rogers  [19] 
previously  reported  a  factor  of  a  hundred  reduction  in  present  in  the  PART  dump  hies 
for  a  previous  version  of  ICEPIC.  This  previous  discrepancy  was  associated  with  the 
two  dimensional  representation  including  a  factor  of  100  difference  in  scale  size  of  the 
simulations.  It  has  been  hypothesized  that  the  current  version  of  ICEPIC  has  had 
this  issue  corrected,  but  a  relic  in  the  old  version  of  the  code  leads  to  the  default 
‘ice.stat’  hie  still  containing  the  factor  of  a  hundred  difference.  This  is  supported  in 
that  the  factor  has  does  not  seem  to  have  an  effect  on  the  modeling  process.  Further 
confirmation  of  this  explanation  for  the  discrepancy  comes  from  the  fact  that  when 
three  dimensions  are  simulated  the  ‘ice.stat’  particle  charge  reported  is  the  same  value 
predicted  by  Equations  83  and  84. 

3.3  Results  and  Analysis 

One  of  the  main  worries  with  window  failure  in  HPM  systems  come  from  the  power 
deposited  on  the  the  dielectric  window.  However,  previous  work  by  Lau  indicates 
the  electrons  are  not  the  primary  cause  of  window  failure,  but  trigger  outgassing 
of  impurities  contained  in  the  dielectric  window  [13].  Ions  would  then  also  become 
a  constituent  of  the  multipactor  plasma  delivering  large  amounts  of  power  to  the 
dielectric  window.  Neuber  and  Edminston  have  published  results  consistent  with  this 
theory  in  their  attempts  to  mitigate  dielectric  window  breakdown  by  preconditioning 
the  dielectric  material  used  in  the  window  prior  to  experimentation  [17]. 

Reflection  by  the  multipactor  electrons  presents  another  problem.  Even  though 
this  was  briefly  studied  in  Rogers  work[19],  reflected  powers  were  taken  over  a  much 
larger  range  of  frequencies,  in  this  research,  giving  a  much  more  complete  param¬ 
eter  space  for  assessing  the  reflection  associated  with  multipactor.  The  results  are 
indicated  in  Table  5. 
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Table  5.  The  maximum  value  of  reflected  power  values  measured  by  the  “probe  line” 
associated  with  Enef.  The  simulation  results  are  compared  to  predicted  values  for 
reflection  calculated  from  equations  73  and  75. 


Frequency  1  GHz 

E(MV/m) 

Observed  Power  Reflection  % 

Theoretical  Power  Reflection  % 

1 

0.04% 

2.04  % 

2 

1.09% 

2.56  % 

3 

4.00% 

2.96  % 

Frequency  2A5GHz 

E(MV/m) 

Observed  Power  Reflection  % 

Theoretical  Power  Reflection  % 

1 

0% 

1.11  % 

2 

0.02  % 

1.81  % 

3 

0.08  % 

2.25  % 

Frequency  10  GHz 

E{MV/m ) 

Observed  Power  Reflection  % 

Theoretical  Power  Reflection  % 

1 

0  % 

0.03  % 

2 

0  % 

0.51  % 

3 

0  % 

0.83  % 

As  indicated  by  the  theory  in  Section  3.1,  the  reflection  of  power  from  the  mul- 
tipactor  is  a  function  of  particle  density,  particle  charge,  the  frequency  of  the  prop¬ 
agating  EM  wave,  and  the  particle  emission  velocity.  As  one  can  see  the  theoretical 
power  reflected  differs  from  the  reported  reflected  power  for  a  few  of  the  simulations. 
The  cause  of  these  differences  may  be  attributed  to  the  particle  distribution  present 
within  the  simulation  being  different  from  the  particle  distribution  associated  with 
previous  distributions  reported  by  Valfclls  [22],  The  theoretical  values  are  based  on 
the  total  count  of  particles  present  within  the  simulation,  which  may  differ  between 
the  two  dimensional  and  three  dimensional  simulation  spaces. 
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IV.  Conclusions  and  Future  Work 


4.1  Conclusions 

The  research  conducted  focused  on  mitigating  the  multipactor  phenomenon  and 
quantifying  the  effects  associated  with  the  reflection  caused  by  the  plasma.  In  or¬ 
der  to  mitigate  the  multipactor,  susceptibility  diagrams  characterizing  multipactor 
were  reproduced  within  ICEPIC.  During  this  process  an  important  discrepancy  was 
observed  for  the  boundary  curves.  The  origin  of  the  discrepancy  was  traced  to  the 
original  method  of  simulation  not  taking  into  account  the  magnetic  field  component  of 
electromagnetic  radiation  propagating  through  the  waveguide  dielectric.  For  the  high 
energy  particles,  the  magnetic  field  imparts  sufficient  momentum  to  return  particles 
to  the  dielectric  surface  even  in  the  total  absence  of  a  bias  electric  field. 

The  remapping  of  the  susceptibility  curves  taking  into  account  the  magnetic  field 
indicates  an  opportunity  to  mitigate  the  initiation  of  multipactor  with  the  applica¬ 
tion  of  a  negative  bias  electric  field.  The  magnitude  of  negative  bias  applied  to  the 
system  would  be  much  less  than  the  corresponding  positive  bias  based  upon  the  lower 
boundary  curve. 

In  quantifying  the  reflection  associated  with  multipactor,  a  new  analytic  treatment 
was  developed  to  determine  the  reflection  coefficient  from  the  particle  density  within 
the  multipactor  plasma.  This  analytic  theory  indicates  that  the  reflection  associated 
with  multipactor  would  be  negligible  in  agreement  with  previous  results  obtained  by 
Fichtl  and  Rogers  [8]  [19].  The  simulation  in  ICEPIC  agreed  for  lower  electric  field 
amplitudes  associated  with  the  plane  wave  normal  to  the  dielectric.  However,  the 
analytic  and  simulated  results  differed  warranting  further  investigation. 
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4.2  Future  Work 


The  research  warrants  further  investigation  in  using  a  negative  bias  as  a  means 
of  multipactor  mitigation.  Further  simulations  need  to  be  conducted  characterizing 
the  new  susceptibility  diagrams  for  various  materials,  taking  into  account  secondary 
emission  properties  associated  with  different  materials.  These  new  susceptibility  dia¬ 
grams  need  to  include  boundary  curves  associated  with  the  lower  and  upper  crossover 
energies  present  in  the  secondary  electron  emission  curve.  New  research  needs  to  also 
investigate  methods  of  implementing  a  negative  bias  field  into  a  real  world  appli¬ 
cations.  The  implementation  needs  to  provide  a  bias  at  a  level  that  will  prevent 
multipactor  initiation  while  also  not  interfering  with  the  propagation  of  electromag¬ 
netic  radiation  through  the  microwave  dielectric.  Spaulding  is  currently  pursuing  this 
study  [21]. 

In  order  to  better  characterize  reflection  in  actual  microwave  systems,  simulations 
also  need  to  be  run  for  time  periods  consistent  with  those  present  in  an  actual  HPM 
system.  Modeling  done  was  usually  performed  over  10  ns  whereas  in  an  actual  system 
microwave  pulses  last  in  the  100-300  ns  scale.  Within  the  simulations  that  considered 
space  charge,  the  total  particle  count  of  the  simulation  was  observed  to  grow  through¬ 
out  the  entirety  of  the  simulation,  indicating  that  power  deposited  on  the  window  and 
reflected  power  would  be  greater  than  values  at  earlier  times.  Characterization  of  the 
overall  particle  growth  as  a  function  of  the  amplitude  and  frequency  of  propagating 
microwave  radiation  would  allow  for  the  forecasting  of  these  values  while  avoiding  the 
computational  requirements  in  simulating  such  large  time  scales. 

Other  mitigation  attempts  should  also  be  explored,  such  as  attempts  at  modifying 
the  geometry  of  the  waveguide  dielectric.  Previous  work  accomplished  by  Rogers  [19] 
in  creating  susceptibility  curves  in  MATLAB®and  Spaulding  [21]  have  indicated  that 
such  attempts  would  prove  fruitful  in  preventing  the  initiation  of  multipactor.  In 
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particular  creating  simulations  within  ICEPIC  characterizing  the  growth  and  decay 
of  multipactor  when  the  incident  plane  wave  is  off  normal  to  the  dielectric  for  the 
space  charge  case  would  be  very  effective  in  showing  the  benefits  of  such  a  mitigation 
method. 

Within  the  research  the  “bowtie”  curves  have  proven  to  be  an  effective  means 
at  understanding  the  time  evolution  of  multipactor  along  the  waveguide  dielectric. 
Looking  at  the  variation  in  the  “bowtie”  curve  for  higher  microwave  field  amplitudes 
to  analyze  the  coupling  effect  with  space  charge  and  reflection  would  also  prove  useful 
in  understanding  the  dynamic  evolution  of  multipactor  in  the  system.  By  observing 
the  electric  field  perpendicular  and  parallel  within  various  areas  of  the  electron  cloud, 
the  spatial  distribution  of  the  reflection  into  the  multipactor  could  also  be  better 
quantified. 
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